Nature of the low temperature ordering of Pr in PrBa 2 Cu 3 6+ 
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Theoretical model is presented to describe the anomalous ordered phase of Pr ions in PrBa2Cu30e+a: below 
Tp r ~ 12 — 17 K. The model considers the Pr multipole degrees of freedom and coupling between the Cu and 
Pr subsystems. We identify the symmetry allowed coupling of Cu and Pr ions and conclude that only an ab- 
plane Pr dipole ordering can explain the Cu spin rotation observed at Tp r by neutron diffraction by Boothroyd 
et al. [A. T. Boothroyd et ah, Phys. Rev. Lett. 78, 130 (1997)]. A substantial enhancement of the Pr ordering 
temperature is shown to arise from the Cu-Pr coupling which is the key for the anomalous magnetic behavior in 
PrBa 2 Cu 3 06+„. 

PACS numbers: 



I. INTRODUCTION 

Non-superconducting Prf^CuaOe+z compound attracted 
considerable attention in the last two decades because of 
its intriguing magnetic and electronic properties within the 
ReBa2Cu306+a; family of compounds (Re=rare earth atom). 1 
In spite of the great efforts from both experimental and theo- 
retical sides, numerous problems remain unresolved such as 
the suppression of superconductivity and the nature of the 
long-range ordered state of Pr sublattice with a unique, about 
an order of magnitude larger than for other rare earth atoms, 
ordering temperature. Although, there exists a discussion 
whether the ground state of PrBa2Cu306+; C is really a non- 
superconducting, insulating material with a magnetically or- 
dered Cu and Pr sublattices, herein we consider this modi- 
fication of PrBa2Cu306+a; as "canonical" and treat it solely 
herein. We note that superconductivity in PrBa2Cu306+ :r was 
reported 2 - and a non-magnetic Pr ground state was found by 
141 Pr NMR.- It was suggested that these observations may 
result from a non-stoichiometric compound, i.e. when Pr oc- 
cupies only about a half of the rare-earth sites, and the other 
half is occupied by the non-magnetic B&£ 

The theoretical model of Fehrenbacher and Rice, i.e. the 
hybridization of Pr 4/ and the nearest neighbor O 2p or- 
bitals is the most accepted model to account for the absence 
of superconductivity in PrBa2Cu3(D6+ :Z ;- It is suggested that 
the localization of holes in the hybridized 4f-2p orbitals ren- 
ders the material non-superconducting. 5 Although the non- 
superconducting nature of PrBa2Cu3C>7 compound is a chal- 
lenging and interesting problem, we concentrate herein exclu- 
sively on the low-temperature Pr ordering. 

PrBa2Cu306+a; is an insulator for every x and the Cu spins 
in the CuC>2 planes order antiferromagnetically at tempera- 
tures of T N - 350 K and 250 K for x = and x = 1, 
respectively^— As the temperature is further decreased, the 
Pr sublattice also undergoes a phase transition?' 10 in the tem- 
perature range Tp r ~ 12 — 17 K depending on x, which ap- 
pears as an anomaly in the temperature dependence of thermo- 
dynamic quantities Jiii 2 - Mossbauer spectroscopy, L1 neutron 
diffraction , 8 ' 11 ' 12 and NMR^i 4 - showed that the Pr magnetic 
moments order antiferromagnetically below Tp r . However, 
there is no consensus among these experiments with respect 



to the magnitude of the ordered Pr moment, its direction, and 
the nature of this transition. 

Neutron diffraction studies found 0.56/^s (x = 0.92) and 
l.lbfiB (x = 0.35) ordered moment with direction tilted 
out from the afo-planeJ. The rotation of Cu moments within 
the CuC>2 plane was found to accompany the Pr ordering* 8 *^ 
which was suggested to result from the strong coupling be- 
tween the Pr and Cu subsystems. However, the magnitude of 
the coupling, its character, and its influence on the magnitude 
of the ordered Pr magnetic moment, and for the ordering tem- 
perature is yet unexplained. 

We give herein a comprehensive description of the Pr order- 
ing in PrBa2Cu30g on the basis of a model of localized 4/ 2 
electrons of the Pr ion. We construct the general form of the 
Cu-Pr interaction allowed by the symmetry and identify the 
Pr order parameter from the rotation of Cu spins below Tp r 
observed by neutron diffraction^^ Although pseudo-dipole 
Cu-Pr interaction has been already proposed j 8 - 15 ' 16 we derive 
the Cu-Pr interaction in general including all Pr multipole 
moments and arbitrary wave vectors. We describe a crys- 
talline electric field (CEF) model which takes a ground state 
quasi-triplet composed by a doublet and a singlet states and 
also includes coupling between the Pr and Cu subsystems. 
This low-energy scheme is consistent with the neutron diffrac- 
tion results 1 ^ and also with the high-temperature susceptibil- 
ity datajii We also discuss the enhancement of the Pr ordering 
temperature due to the Cu-Pr coupling and the magnitude of 
Pr magnetic moment in the ordered state. 



II. GENERAL FORM OF THE CU-PR INTERACTION 

In PrBa2Cu3(D6, the Cu spins are antiferromagnetically or- 
dered along the [100] direction in the Q1O2 planes with the 
ordering vector Q = [ir/a, tv/a, ir/c] below the temperature 
Tcu ~ 350 K (phase AFI)4 As the temperature is further de- 
creased, the Pr sublattice also undergoes an ordering at tem- 
perature Tp r = 12 K. This ordering is accompanied by the ro- 
tation of Cu spins within the CuC>2 planes in such a way that 
the new spin structure is non-collinear along the c-direction 
(phase AFIII)i 8 ' 15 ' 18 The Cu spin rotation is characterized by 
the ordering vector Q = [w/a, 7r/a, 0]. The Cu spin structure 
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FIG. 1: Magnetic structure of Cu spins in the Cu02 planes and Pr 
dipole moments in the ab-plane in the temperature range T < Tp r . 
The central Pr ion in the a6-plane is located at position ro = 
(0,0,0), and the numbers label the Cu sites, i.e. Cu[l-8] used in 
the text. 

in the C11O2 planes is shown in Fig.[T]in the temperature range 
T < T Pr . 

First, we construct the invariant form of the interaction be- 
tween the multipole moments of a Pr ion at position ro = 
(0, 0, 0) and the surrounding eight Cu spins, Cu[l-8], shown 
in Fig.Q] Although the local symmetry at Pr site is tetragonal, 
a Cu-Pr pair has lower symmetry. For example, the Cu[l]-Pr 
pair has reflection symmetry with respect to the [1,1,0] mir- 
ror plane. The transformation of total angular momentum J 
of Pr ion under this operation is: J x — > J y , J y — > J x , and 



J z — > — J z . The same transformation is applied for the Cu 
spin S. Thus, spin operators J x — J y {S x — S y ) and J z (S z ) 
are odd, while J x + J y (S x + S y ) is even under this reflection. 
Not only the bilinear products of the Pr dipole moments and 
Cu spins appear in the interaction but also the bilinear prod- 
ucts of rank-3 octupole and rank-5 triakontadipolei2 operators 
of Pr ion and Cu spins since they are allowed by the time- 
reversal symmetry. The octupole operators Tf and T@ — T,^ 
are even, while T xyz and T% + T@ are odd operators under 
the present transformation. In the case of tetragonal symme- 
try, there is one more independent magnetic operator, namely 
the rank-5 triakontadipole operator V\ u = J x J y J z (J x — J y ) 
which is even under the present reflexion. Thus, the invariant 
form of the interaction for the Cu[l]-Pr pair is constructed as 

^? u " Pr [l] = (S x + S y ) [ Cll {J x + J y ) + c 12 T? + c 13 A lu ] 

+ (S x - Sy) [C2l(j x - Jy) + C2%J Z + C 23 T xyz ] 
+ S z [c 31 (J x - J y ) + C 32 J Z + C3 3 T xyz ] , (1) 

where cy are constants which are not determined by the sym- 
metry itself. The pair Cu[2]-Pr is connected to the pair Cu[l]- 
Pr by a 7r/2 rotation around the c-axis. Under this operation 
the transformation of the angular momentum components is 
given by J x — > — J y , J y — > J x , and J z — > J z . The same 
transformation holds also for the Cu spin components. Thus, 
the form of the interaction for the Cu[2]-Pr pair can be ob- 
tained from Eq. ([T]i by applying the ir/2 rotation. Repeating 
further the appropriate rotations and reflections, the interac- 
tion can be obtained for all the eight Cu-Pr pairs. The Cu 
spin at position r is expressed by Fourier transformation as 
S(r) = J2q S q c' qr . Using this expression, we obtain the to- 
tal interaction as 



H Cu-Pr = ^ H Cu-Pr [fc] = 8 J] {(c n + C 21 )(S^J X + Sp y ) Cx C y C z + ( Cll ~ C 21 )(S^J y + S*J x )s xSy C z 
k=l q 
+ c 22 J z [S^s X CyS z SyC X SyS Z ^ + C31 [J X S^S X SyC Z JyS^C X SyS Z } + c 32 J z S* C X Cy C Z 

S X CyS z "J" O C X SyS z ] + c 13 A lu [S% Sz + S^SxCySz] + c 23 T xyz [S^c X SyS Z S^S x CyS z ^ 

+ c 33 T xyz S^s x s y c z } , (2) 

I 



where Cfc(sfc) denotes cos(q , fea/2)(sin((/fca/2)) for k=x, y, 
and cos(<7fcc/2)(sin(<7fcc/2)) for k=z. 

Now we discuss the multipole order of Pr sublattice based 
on expression @. In phase AFI, the Cu spin component S* q 
with q = Q is non-zero. In this phase there is no cou- 
pling between the Cu spins and Pr magnetic moments be- 
cause the wave vector Q leads to c x = c y = c z = and 
s x = s y = s z = 1, which gives that none of the terms sur- 
vives in expression (O. On the other hand, a coupling of Cu 
spins and Pr moments is present in phase AFIII. The non- 
collinear Cu spin structure in phase AFIII can be described 
by the appearance of the extra spin component in addition 



to the component The wave vector Q gives c, x = c y = 0, 
c z =l, and s x = s v = 1, s z = in expression ©. Together 

with the condition 7^ 0, we find that the only surviving 

term in the form of Cu-Pr interaction is the coupling J x . 

We note that the coupling S^J y is also allowed by symmetry 
since the [100] and [010] domains in phase AFI are equiva- 
lent, which gives that the domains (S?) ^ 0, (S^) = and 
(S$) + 0, (S$) = are also equivalent. 

Thus, we conclude that the neutron diffraction result for 
the Cu spin rotation below Tp r is compatible with Pr dipole 
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moments lying along the crystalline [100] or [010] direction, 
and no deviation from this direction is allowed by symmetry in 
contrast to the structure proposed in Ref.[|l Since the (S$) 
Cu spin components alternate along the a and b directions, 
the Pr dipole moments also form an AFM structure in the ab- 
planes with the ordering vector Q. The magnetic structure 
for the Cu spins and Pr dipole moments below Tp r has the 
structure shown in Fig.Q] 24 



III. CEF MODEL FOR THE COUPLED CU-PR SYSTEM 

We now turn to discuss a CEF model of the Pr ion 4/ 
electrons to describe the Pr ordering in PrE^CuaOg. The 
model consists of the coupled Pr and Cu subsystems and al- 
lows to calculate experimentally relevant quantities. The Pr 
ion has 4/ 2 electronic configuration in PrBa2Cu30g, which 
gives J = 4 as the total angular momentum based on the 
Hund's rule for the ground state. In the Pr subsystem, we con- 
sider AFM interaction between the in-plane J x and J y dipole 
moments described by the Hamiltonian 



H 



Pr 



— 2 A [J x ,iJ x ,j + Jy,iJy,j] 



(3) 



(hi) 



where z = 4 is the number of nearest-neighbor Pr ions and 
A is the coupling constant. The Hamiltonian <£3j is taken in a 
quasi-triplet subspace, where the ground state is the doublet 



\d±) =a|±3) + v / l-a 2 |Tl> 
with a near singlet excited state 



\s) 



V2 



2)), 



(4) 



(5) 



where the states are expressed by the eigenstates of J. The 
doublet and singlet states are separated by the gap AM*H The 
other CEF levels lie at much higher energy therefore these can 
be neglected when describing the low temperature behavior. 

We treat the Cu subsystem phenomenologically and con- 
sider the free energy expansion: 



-^Cu = --lnTr s (c 



= a s ((S?) 2 + is?) 2 ) + b s ((S?) 2 + (£^) 2 ) 2 
+ 4 Cs (S«) 2 (S,p) 2 +J-° u , (6) 

where J 7 ^ contains the part corresponding to the Cu spin 
components, and S~*, which are not relevant for our 

discussion. We assume that (S^LO is not critical around 
T = Tp T , thus, a s > 0. Parameters b s and c s determine 
the anisotropy as we discuss below. 

The local coupling between the Cu and Pr subsystems is 
assumed as 



H c = z'Y,x[j IB ,i(S?) + JvAS2) 



(7) 



which corresponds to the only non-vanishing term in expres- 
sion (f2|i for Q = [it /a, ir/a, 0]. Here, z' = 8 is the number 
of Pr nearest-neighbor Cu ions and A is the Cu-Pr coupling 
constant. The series expansion of the total free energy has the 
form 



T = 



where 



-ilnTr s 



-mcu 



) - ^Tr, (c 



- / 3(w Pr +-H c ; 



Cm 



(8) 



Jo = a J ((J ; T r > 2 + (./ y Pr ) 2 ) + 6 ./((^ Pr > 2 

+ (J^-) 2 ) 2 + \z'({J^){S$) + (J^)(S?)) 
+ a s ((S?) 2 + (S?) 2 ) + b s ((S?) 2 + (fij) 2 ) 2 
+ 4c s (5^) 2 (5«) 2 . (9) 
The coefficient aj is expressed in Eq. (O as 



_ 1 2 /l 1 (2 7 ) 2 (e-^-l) 
° J 2 U A (e-^ + 2) 



(10) 



where 7 = 1/4(9 — 2a 2 + 6\/7a\/l — a 2 ) is the matrix ele- 
ment of the dipole operator, J x , between the doublet and sin- 
glet states. 

The quasi-triplet subspace does not carry anisotropy with 
respect to the dipole moment components J x and J y , thus the 
Pr magnetic moment is isotropic in the a6-plane for the Pr 
subsystem. However, there is a weak anisotropy even for the 
structurally tetragonal x = compound which gives rise to a 
degenerate [100] and [010] easy axis ordering of the Cu spins 
in the AFI phased We introduce this anisotropy of the Cu 
sublattice through the parameters b s and c s in expression (O, 
which produces anisotropy also for the in-plane Pr magnetic 
moments through the Cu-Pr coupling. In the AFIII phase, 
these easy axes are inherited and the Cu spin ordering with 
wave vector Q is reproduced when c s > 0. We assume an in- 
finitesimally small magnetic field along the [010] direction on 
the Cu sublattice, which resolves the [100], [010] degeneracy 
and stabilizes the solution (J x ) 7^ 0, ( J y ) = against the so- 
lution ( J x ) = 0, (J y ) 0. The resulting Pr dipole moments 
have also a [100] easy axis. 

In our scenario, the Pr dipole moments order at T = Tp T , 
and the order (J x ) ^ induces the (S^) Cu spin component. 
From the condition dJ /d(S^) =0we obtain 



(11) 



Substituting this expression into the free energy expansion T, 
the second order coefficient of ( Jl >r ) is obtained as 



aj = aj 



4a s 



(12) 



We define the temperature To which corresponds to the 
non-coupled system with A = and the temperature Tp r for 
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FIG. 2: Enhancement of the transition temperature, Tp r /Tb, due to 
the coupling parameter A. Inset shows the Cu spin-rotation angle <j> 
at T = as a function of A. The used parameter values are A = 4 K, 
A = 1.2 K, a = 0.96, and a s = 140. Arrows show the experimental 
values for cf> and A from Refs.^and [TJ, respectively. 

the coupled system A ^ 0. To is obtained from the condi- 
tion a j = 0, where the second-order coefficient has the form 
a j = a'j(T — Tq) in the vicinity of the phase transition. How- 
ever, this transition temperature is modified due to the cou- 
pling to the Cu subsystem (A ^ 0). Tp r is obtained from the 
condition aj = 0, which gives 

T Pr =T + ^-. (13) 

Thus the Pr transition temperature is enhanced due to the Cu- 
Pr coupling irrespective of the sign of the coupling parameter 
A. 

The rotation of the Cu spins in the CUO2 plane is expressed 
as <j> = tan _1 ((5^)/(5 1 ^)), where <f> = corresponds to 
the x-direction. The Cu magnetic moment is expressed as 

MCu = .9CuM(SQ) 2 + (S?> 2 ] 1/2 , where g Cu « 2 is the 
Cu g-factor. An ordered magnetic moment of 0.66/WCu[2] 
is found also for the AFI phase due to quantum fluctuations. 21 
Our model includes this effect phenomenologically as this 
ordered magnetic moment is retained above Tp r in the AFI 
phase. 

Figure [2] shows the enhancement of the Pr transition tem- 
perature due to the Cu-Pr coupling and the corresponding Cu 
spin rotation angle as a function of the coupling parameter 
A, calculated within the above model. We used A = 4K 
and a = 0.96 for the CEF parameters, where the latter value 
is obtained from the high-temperature susceptibility datai^ in 
PrBa2Cu3C<6 . Neutron scattering found the interaction param- 
eters as A = 1.2 K and A = 1.74 Ki^ and an ordered Cu mo- 
ment of ncu = 0.64/j,b for x = 0.35, 8 which were used in 
the calculation of Fig.|2]as fix parameters. Clearly, the model 
accurately reproduces the <j> = 20° rotation of the Cu spins 
which was observed experimentally in Ref.[H In addition, the 
model accounts for about 70 % enhancement of the Pr order- 
ing temperature with this parameter set. 
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FIG. 3: Transition temperatures Tq and Tp r as a function of A. The 
parameter values are chosen as A = 1.2 K, A = 1.74 K, a — 0.96. 
Solid line corresponds to second-order, while dashed line to first- 
order phase transition. The ratio Tp r /Tb is also shown as continuous 
line. 



The ground state Pr magnetic moment and the enhance- 
ment of Tp r strongly depends on the value of the CEF en- 
ergy gap, A. In Fig.[3j we show the calculated Tp r as a func- 
tion of A, where we fix the interaction and CEF parameters 
as A = 1.2K, A = 1.74K, a = 0.96. The value of the 
phenomenological parameter a s is chosen for each A value 
so that it reproduces the experimental values <j> = 20° and 
MCu = 0.64^b- For small values of A the phase transition is 
second-order. The transition temperatures To an d Tp x are sup- 
pressed with increasing energy gap A, and the phase transition 
changes to first-order above a critical value of A. Increas- 
ing further A, the phase transition disappears. These behav- 
iors are due to the fact that the J x dipole order is interaction- 
induced within the quasi-triplet subspace since the dipole op- 
erator J x has non-vanishing matrix element only between the 
doublet and singlet states. Around A ~ 10 K the ratio Tp t /Tq 
is considerably enhanced due to the Cu-Pr coupling. 



IV. MEASUREMENT OF THE LOCAL FIELDS 

Finally, we discuss the field-angle dependence of the local 
magnetic field acting on a Pr ion. This local field could be 
measured by a local magnetic probe such as e.g. 89 Y using 
NMR 2 * or Gd 3+ using ESRS which can be substituted into 
the Pr sublattice in a low concentration. 

The magnetic field H pro b c acting on the local magnetic 
probe has five different sources: exchange and dipole fields 
from the surrounding Cu spins, exchange and dipole fields 
from the surrounding Pr dipole moments, and the external 
magnetic field. Contributions except the external magnetic 
field are directed along the ^-direction, and we express their 
effect by the field h x . We keep the field h x fixed as the direc- 
tion of the external magnetic field H is changed which is the 
situation for small values of the magnetic fieldi 2 ^ For a gen- 
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eral direction of the external magnetic field H, we write the 
Hamiltonian of the local probe as 



^probe — Tl'yS • Hp ro bc 



(14) 



where 7 is the gyromagnetic ratio of the local probe, and S 
is its spin (either electron or nuclear). Thus, we obtain the 
experimentally detected magnetic resonance shift, /i pro bo (in 
magnetic field units): 



hprobc — S(hj- os ho), 



(15) 



where h re s is the resonance field and ho is the reso- 
nance position for the AFI phase. The external mag- 
netic field is expressed by field angles as H/|H| = 

[H x ,H y ,H z ] = [cos0sin0,sin</>sin0, cos0], and we as- 
sume that [S x , S y , S z ] ~ [H x , H y , H z ]. This gives 



h pI 



probe 



{h 3 
H(cos6Y 



H cos d> sin 9) cos d> sin ( 



H(s'md>sm( 



(16) 



Figure|4]shows the field angle dependence of /i pro be by chang- 
ing the direction of the magnetic field in the planes [001] and 
[110]. We estimate the change of /i pro bc as 



"-probe 

[100]- Vobdiio] 



(17) 



frprobc[110] - /lprobo[001] 

which ratio can be checked in the NMR or ESR experiments. 

V. DISCUSSION AND SUMMARY 

We studied the nature of the Pr ordered phase in 
PrBa 2 Cu 3 06 in a quasi-triplet CEF model of 4/ 2 electrons, 
where we also included the symmetry allowed coupling be- 
tween the Pr and Cu ions. The reason to include Cu-Pr cou- 
pling is the Cu spin rotation observed at the Pr ordering tem- 
perature in neutron diffraction, 815 which has not been ob- 
served for other rare earth ReB^CusOe compounds. In the 
4/ rare-earth series, the localized electrons of Pr ion tend to 
hybridize rather strongly with conduction electrons as in e.g. 
PrFe4Pi2. This can be the reason why the Cu-Re coupling 
is much stronger in ReB^CusOg for Re=Pr compared to the 
other rare-earth ions, Re. 

We showed by a general symmetry analysis that there is no 
coupling between the Cu and Pr subsystems in phase AFI, but 
a coupling emerges in phase AFIII. We prove that only the 
AFM ordering of J x (J y ) dipole moments of Pr ions is consis- 
tent with the Cu spin rotation observed at T = Tp r by neu- 
tron diffraction^^ The interpretation of the neutron diffrac- 
tion data in Ref.H includes a tilting of the ordered Pr magnetic 
moments out of the a6-plane. However, the presence of non- 
zero J z dipole component does not follow from the general 
symmetry analysis as it was shown in Sectionllll 

Upon identifying of the Pr order parameter, we studied the 
Pr ordering temperature in the Cu-Pr coupled quasi-triplet 
CEF model by changing the interaction and CEF parameters. 
We found that the Cu-Pr coupling enhances the Pr ordering 
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FIG. 4: Field-angle dependence of the local field h plo b e acting on 
the local magnetic probe by changing the direction of the external 
magnetic field in the [001] (top) and [110] (bottom) planes with 6 = 
7r/2 and cj> = 7r/4, respectively. In the plot we take the exemplifying 
values h x — 1 and H = 0.5. The field angles are defined in the text. 



temperature, Tp r , compared to the uncoupled ordering tem- 
perature, Tq. 

By fixing the Pr-Pr and Pr-Cu interaction parameters to the 
values obtained by neutron diffraction, 15 and the CEF pa- 
rameter to that obtained by high-temperature susceptibility 
measurements^ we found that the Pr ordering temperature 
is considerably enhanced due to the Cu-Pr coupling in the pa- 
rameter range A ~ 8 — 10 K. Namely, the enhancement of the 
ordering temperature is as large as Tp r /Tn ~ 8 — 12 in this in- 
terval. This observation explains the uniquely large ordering 
temperature for Re=Pr in the series ReBa2Cu306+ 2 > 

In addition, we predicted the magnetic field angle depen- 
dence of the local magnetic field acting at a Pr site in the AFM 
ordered phase of the Pr ions. This can be directly compared to 
measured NMR or ESR spectra which detect the local mag- 
netic field by a local magnetic probe. 

Our model contains the minimal number of parameters 
which is required to account quantitatively for the experimen- 
tal data. Given that the parameter values are assumed from 
independent measurements and no fit is performed, the agree- 
ment is reasonable. First, the realistic estimate A ~ 8 — 10 K 
gives a Tp r of 5 — 3.6 K which is to be compared to the exper- 
imental Tp r of 12 K. We note that this estimated range gives 
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the Pr transition temperature To for the uncoupled system as 
0.7 — 0.3 K, which is close to the typical values of the order- 
ing temperatures for other rare-earth ions such as Yb (0.35 K), 
Nd (0.5 K) or Dy (IK). This fact may also indicate that the 
uniquely large ordering temperature for Pr arises from the Cu- 
Pr coupling. Second, the interval A ~ 8 — 10 K is not very far 
from the estimate A « 17 K based on the analysis of the high- 
temperature susceptibility data. 17 Furthermore, the estimated 
interval for A and also the splitting pattern of the quasi-triplet 
subspace is consistent with the results of neutron scattering, ' 1 
which found that the splitting of the ground state quasi-triplet 
does not exceed 2 meV above Tp r , and two inelastic lines are 
observed at energies 1.7meV and 3.4meV at T = 5K in the 
ordered phase. Finally, for the interval A ~ 8 — 10 K our CEF 
model gives the ordered Pr magnetic moment at low temper- 
atures as jUp r ~ 1.8 — 1.7/iB ) which is to be compared to the 
observed value /ip r = 1.15/Ub.— The quantitative discrepan- 
cies between the results of our model and the experimental 
data arise because of the simplicity of our model and minimal 



number of parameters. 

In summary, the model we propose herein describes the 
main features of PrBa2Cu3C>6 system, and demonstrates the 
importance of the Cu-Pr coupling which may be the key point 
to understand the anomalous behavior of this compound such 
as the enhancement of the Pr ordering temperature. Further 
studies which include also the itinerant character of 4/ 2 elec- 
trons are necessary to give an extensive description of this 
rather complicated system. 
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